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Role of Spin-Orbit Coupling on the Spin Triplet Pairing in NaxCo02 ■ 2/H2O 
II: Multiple phase diagram under the Magnetic Field 
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The possibility of multiple phase diagram in the novel superconductor NaxCo02 • J/H20 
is analyzed on the basis of the multi-orbital Hubbard model including the atomic spin-orbit 
coupling. We have shown that the spin triplet pairing state is stable in this model. The p-wave 
(/-wave) state is stabilized when the Hund's rule coupling is large (small). In the precedent 
paper, we have determined the direction of d-vector at T = Tc and H = within the linearized 
Dyson-Gorkov equation. In this paper, the pairing state below Tc and under the magnetic field is 
determined within the weak coupling approximation including the paramagnetic effect. We find 
that the p+f coexistent state is stabilized at low temperatures in a part of parameter range. We 
point out that the phase diagram in the H-T plane is quite different between the p-wave, /-wave 
and p + /-wave superconductivities. The characteristics of each phase are clarified by showing 
the magnetic susceptibility and specific heat. We discuss the comparison with experimental 
results and suggest some future experiments to detect the multiple phase transition. 
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1. Introduction 

The spin triplet superconductivity is one of the most 
exciting topics in condensed matter physics. Vast stud- 
ies have been devoted to the candidate materials such 
as Sr2Ru04,^) (TMTSF)2PF6,2) UPta,^) UNiaAla,*) 
UGe2^-' and URhGe.®^ Recently, a novel superconduc- 
tivity has been discovered in NaxCo02 • yH20 ''^ and the 
possibility of spin triplet superconductivity has attracted 
huge interests. 

The properties of superconducting state in NaxCo02 ■ 
yH20 have been investigated by many experimen- 
tal studies including the mag netization,^) NMR,^"^^) 
/iSR,^"*"^®^ specific heat,^^"^^^ resistivity,^*^^ and impu- 
rity effect. Except for the impurity effect, most of 
them indicate a non-s-wave pairing. The discovery of 
magnetic phase in the neighborhood of superconducting 
phase^^' ^'^^ indicates a strong electron correlation which 
generally leads to an unconventional superconductivity. 

As summarized in Ref. 24, there have been many theo- 
retical studies. Although most of these studies have 
assumed single-orbital models, the justification of this 
assumption is not clear because the conduction band 
in NaxCo02 • yH20 has orbital degeneracy, as pointed 
out by Koshibae et al.^^^ In order to examine the super- 
conductivity in the multi-orbital system, we have con- 
structed a three-orbital Hubbard modeP^^ which appro- 
priately reproduces the electronic structure obtained in 
the LDA calculations.^^"'"'-' From the results of pertur- 
bation thcory^^^ and FLEX approximation, '^^^ we have 
shown that the spin triplet superconductivity is stable 
in the wide parameter range. Then, the p-wave super- 
conductivity and /-wave superconductivity are nearly 
degenerate because of the character of orbital in each 
Fermi surface.'*^-' We have found an orbital dependent 
superconductivity like Sr2Ru04^^^ and derived the two- 
orbital Hubbard model reproducing the e' -doublet which 



mainly leads to the superconductivity. 

In the precedent paper, ^''^ which we call "I" in the 
following, we have constructed a two-orbital Hubbard 
model including the spin-orbit coupling term and deter- 
mined the d-vector at T = Tc with use of the linearized 
Dyson-Gorkov equation. In the present paper, we study 
the phase diagram in the temperature and the magnetic 
field plane (H-T plane) and clarify the physical proper- 
ties, such as magnetic susceptibility and specific heat, in 
each phase. This study is for the following two purposes. 
(1) One purpose is to provide clear physical quantities 
which can be compared in experiments. Many measure- 
ments including the NMR Knight shift are performed 
under the magnetic field. Since the d-vector may rotate 
under the magnetic field so as to gain the Zeeman cou- 
pling energy, it is highly desired to understand the d- 
vector under the magnetic field. (2) Another and more 
interesting purpose is to suggest a possibility of multiple 
phase transition in the superconducting state. The mul- 
tiple phase diagram may appear in the H-T plane owing 
to the rotation of d-vector. If the multiple phase tran- 
sition is observed in experiments, that will be a strong 
evidence for the spin triplet pairing. 

The effects of spin-orbit coupling ware clarified in I. 
In the model without any spin-orbit coupling term, the 
d-vector rotates in the infinitesimal magnetic field. In 
general, the d-vector is determined by the competition 
between the magnetic field and the anisotropy arising 
from the spin-orbit coupling. For superconductors, it is 
reasonable to consider the atomic spin-orbit coupling, 
namely so-called L-S coupling as a microscopic origin of 
the spin-orbit coupling for Cooper pairs. Therefore, the 
d-vector in spin triplet superconductors is a fundamental 
issue of multi-orbital systems. 

The multiple phase transition in heavy fermion su- 
perconductors has attracted much interests for last two 
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decades. ^•'""'^^ Then, the theoretical treatment has rehed 
on the phenomenological theory and the role of L-S cou- 
pling on the rf-vector has been a longstanding problem. 
For instance, the pairing symmetry in UPta has been 
discussed for a long time and the anisotropy of d-vector 
is still under debate. We consider that the study 
on the d-electron systems will be a first step for this 
problem. The microscopic theory including the spin-orbit 
coupling can be applied to the d-electron systems like 
NaxCo02 • 2/H20 and Sr2Ru04 because these systems 
have a simple electronic structure in comparison with 
/-electron systems. Actually, we have performed a mi- 
croscopic analysis on the d- vector in Sr2Ru04'' "^^^ and 
obtained the consistent results with experiments.^^' 
We have found that Sr2Ru04 is a particular case in the 
sense that the anisotropy of d- vector is very small owing 
to the orbital dependent superconductivity. This result 
has been confirmed by the NMR measurement. In this 
paper, we show that NaxCo02 • y^2^ is a more general 
case where the spin-orbit coupling plays a more impor- 
tant role. We point out that the multiple phase transition 
will occur when the order parameter includes the p-wave 
component. 

This paper is organized as follows. In §2, we briefly 
summarize the results of linearized Dyson-Gorkov equa- 
tion for the two-orbital Hubbard model and derive an 
effective single-band model including the Zeeman cou- 
pling term. In §3, we determine the phase diagram in H- 
T plane within the mean field theory. We obtain three 
pairing states at zero magnetic field, namely the p-wave, 
/-wave and p -\- /-wave states. The vector under the 
magnetic field are investigated for these states in §3.1, 
§3.2 and §3.3, respectively. Our interests are focused on 
the magnetic field parallel to the two-dimensional plane 
because a wide variety of phase diagram is expected 
in this direction and also because the Knight shift has 
been measured under the parallel magnetic field in most 
cases.^'^°'^'^'^'*'^^'^^) The nature of multiple phase tran- 
sition is analyzed with attention to the thermal and mag- 
netic properties. In §4, the role of vortex state and that 
of aig-orbital which are ignored in §3 will be discussed. 
We show that an extra phase transition may be induced 
by these effects. Some discussions are given in the last 
section §5. 

2. Formulation 

2. 1 Pairing state at T = Tc and H = 

As in I, we adopt a two-orbital Hubbard model repro- 
ducing the eg-doublct. From the analysis of three orbital 
Hubbard model without including the spin-orbit cou- 
pling, we have found that the superconductivity is mainly 
induced by the Cg-orbitals."^^''*^^ Indeed, the two-orbital 
Hubbard model without including the aig-orbital appro- 
priately reproduces the results of three-orbital Hubbard 
model. Wc choose the parameters of two-orbital Hub- 
bard model as in I unless we specify. 

The pairing symmetry in the two-orbital and three- 
orbital Hubbard models are basically determined by the 
two parameters, namely Ue and Jh- Here, tie is the holes 
in the Cg-Fermi surface and Jh is the Hund's rule cou- 
pling. We have found that the /-wave superconductivity 
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is favored when rio is large and vice versa.'^^'^^'^ This is 
because the amplitude of order parameter in the p-wave 
symmetry have to be small around the K-point owing to 
the periodicity of Brillouin zone. As increasing rie, the 
eg-Fermi surface approaches to the K-point and the p- 
wave superconductivity becomes unfavorable. As for Jh, 
the p-wave (/-wave) superconductivity is stabilized when 
Jh is large (small). 

The rf-vector is determined by the spin-orbit coupling 
term which plays an essential role for the issues in this 
paper. The coupling constant 2A of spin-orbit coupling 
term has been estimated as 57meV from NMR measure- 
ments.^"^) This value corresponds to A = 0.17 in our unit 
if we choose the band width W ^ 9 = 1.5eV according 
to the LDA calculations. Although this estimation has 
some ambiguities, the spin-orbit coupling in this order is 
much smaller than the band width. This is in contrast 
with the RVB theory discussed by Khaliullin et al.^^^ 
who have considered the opposite limit A = oo. 

If the spin-orbit coupling is neglected, there remains 
a 3 X 2 = 6-fold degeneracy in the p-wave state, while 
the degeneracy in the /-wave state is 3-fold. Due to the 
spin-orbit coupling, the degeneracy is lifted and the pair- 
ing states are classified into Pxy+, -fxy- and for the 
p-wave symmetry while F^y and for the /-wave sym- 
metry. The d-vector in these states has been summarized 
in Table II of I.^^^ We have determined the pairing state 
at T = Tc and H = with use of the linearized Dyson- 
Gorkov equation whose derivation for the SU (2) asym- 
metric system has been given in I and Ref. 47. For the 
two-orbital Hubbard model, we have obtained the two 
independent equations as, 

AlTAn(fc) = -^yTT(fc,fc')|G2(fc')pATT(A='), (1) 

fe' 

xl-^Anik) = -J2ynik,k')\G2ik')\'An{k'). (2) 
fe' 

Here, A|| and A|| arc the order parameters for the equal 
spin pairing state and opposite spin pairing state, re- 
spectively. For the effective interactions, V^i{k,k') and 
Vi;i{k, k'), we applied the second order perturbation the- 
ory (SOP) with respect to the Goulomb interactions, i.e. 
U, U', Jh and J. Note that we have confirmed that 
the SOP is consistent with the renormalized third order 
perturbation theory as well as the FLEX approximation 
which include higher order corrections.^^' The Green 
function G2{k) is defined as G2{k) — -"53 (fc) ^^^^^^ 
E2{k) is the upper band composing the Fermi surface. 
Note that E2{k) is slightly affected by the spin-orbit cou- 
pling. 

By solving cqs. (1) and (2), we have found that the 
spin-orbit coupling stabilizes the Pxy-i--state when the 
pairing symmetry is p-wave. Then, the pairing state 
is described by the d-vector as, d = p^x + Py-y or 
d = PyX — p^y. There is no violation of time-reversal- 
symmetry in this state which is consistent with the /iSR 
measurement.^'*) In the case of /-wave symmetry, F^y- 
state (d = fix — af2y or d = af2X + fiy) and -F^-state 
{d = fiz) can be stabilized depending on rig and Jh- 

As shown in I, the anisotropy of d-vector defined as 
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|AP — AJ^|/A|^ is much larger in the p-wave state than 
in the /-wave state. This is because the first order term 
in spin-orbit coupUng vanishes in the latter case. This 
means that the d-vector in the /-wave superconductiv- 
ity is rotated by a small magnetic field parallel to the 
d- vector. Thus, the role of spin-orbit coupling is quite 
different between the j;-wavc! and /-wave superconduc- 
tivities. This leads to the qualitatively different phase 
diagram under the magnetic field as shown in §3.1-3. 

2.2 Weak coupling approximation 

In order to determine the pairing state below Tc and 
under the magnetic field, we introduce an cfFectivc single- 
band model which reproduces the results of linearized 
Dyson-Gorkov equation at H = 0. 

k,s 

-^91 01x(fc)</'lx(fc')cI,scLfe,sC-fe',sCfe',^. 
fc,fe',s 

-^92 01y(A;)01y(fe')c^cLfc_^C_fe',sCfc',s 
fc,fe',s 

-^93 Y [i«'/'ix(fe)'?^iy(fe')4,scLfe,sC-fc',«Cfe',s 

k,k' ,s 

-is0iy(A;)?iix(fe')4,scLfe,sC-fe',sCfe',s] 

-^94 Y 02x(fc)</>2x(fe')c^clfc _^C_fc/ _sCfe/,s 

fe,fc',S 

-^95 Y ^2y{k)(p2yik')cl^^c'_^_^C-k',-sCk',s 
k,k' ,s 

-MH, (3) 

where c\. ^ is the creation operator for the ii'2(fe)-band 
electron with pseudospin s. The coupling constants gi 
are determined so as to reproduce the result of linearized 
Dyson-Gorkov equation, namely eqs. (1) and (2). There- 
fore, these parameters are determined by the microscopic 
parameters of Hubbard model such as rig and Jh- Here, 
4>ix{k) and 4>iy{k) arc the orbital part of Cooper pair 
wave function for the d-vector along the plane which 
are determined by the linearized Dyson-Gorkov equa- 
tion in eq. (1). For instance, we obtain A||(fc,i7rT) = 
— 0ix(fe) + i(/)iy(fe) when the Pxy+-state is stabilized at 
T = Tc. Similarly, (^2x(fe) and <p2y{k) are the wave func- 
tion for the (i- vector along z-axis determined by eq. (2). In 
case of the j?-wave superconductivity, <pix{k) and (/'iy(fe) 
(i = 1,2) denote the wave functions in the Px- and py 
wave symmetry, respectively. In case of the /-wave su- 
perconductivity, these are wave functions in the /2- and 
/i-wave symmetry, respectively. We have neglected the 
frequency dependence of (f>ia{k) for simplicity and use 
the value at u;„ = inT, i.e. the smallest Matsubara fre- 
quency. 

As a role of the magnetic field, we take into account 
the Zeeman coupling term, namely the last term in 
eq. (3). This term plays an essential role for the rotation 
of d-vector. We obtain the magnetic moment from the 
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two-orbital Hubbard model as Ma = {La + 2Sa)HB = 
X^/fe 5a(fc)'§a(fc)A*B where Sa is the operator of pseu- 
dospin for the £'2(fc)-band electron. Here, the momen- 
tum dependent ^-factor is obtained as ga{k) = 2 < 
k,Sa = \\La + 2Sa\k,Sa = 5 > whcrc \k,Sa = i > 
is the wave function of £^2(fc)-band with Sa = ^. Note 
that a = X, y and z is the symmetry axis of crystal and 
g^{k) = gy{k) ^ gz(k). In the following, we consider the 
magnetic field along x-axis unless we mention. 

In the following calculation, we simply ignore the or- 
bital effect which induces the vortex state. We discuss the 
orbital effect in §4.1 and show that the phase diagram in 
the high field region H Hc2 can be affected when the 
orbital part of order parameter has multi-component as 
in the p-wave superconductivity. 

Because we have assumed separable pairing interac- 
tions, the following two effects are ignored in eq. (3). 
First, although the momentum dependence of wave func- 
tions (j)ia{k) {i = 1,2 and a = x, y) depends on the 
temperature and magnetic field in general, ''^^ we have 
ignored this deformation of (pia{k). Second, the feedback 
effect, namely the effect of superconductivity on the effec- 
tive interaction Vss'{k, k') is neglected. Although a part 
of the feedback effect can be represented by the tem- 
perature dependence of attractive interaction, we have 
introduced temperature independent coupling ccmstants 
gi in eq. (3). These simplifications are justified in the 
weak coupling region T^/W -C 1 which is reasonable in 
NaxCo02 • J/H2O. Qualitatively the same assumptions 
have been adopted in the analysis of Sr2Ru04.^'') 

In principle, these effects can be taken into account 
in the non-linear Dyson-Gorkov equation for the two- 
orbital Hubbard model including the Zeeman coupling 
term. However, the small value of transition tempera- 
ture Tc ~ 5K makes the numerical calculation very dif- 
ficult. Note that we have to calculate at T ^ 0.5K in 
order to determine the pairing state at T = O.lTc. The 
extraordinary large size simulation is needed to discuss 
such a small energy scale. If a larger value of Tc than 5K 
is assumed to enable the numerical calculation, we will 
overestimate contributions ignored in the weak coupling 
theory. For instance, the renormalization of effective in- 
teraction due to the superconducting order parameter 
is overestimated. As for the effect of magnetic field, the 
contributions scaled by (ibH/W are also overestimated, 
although only the contributions scaled by (ibH/Tc are 
taken into account in the weak coupling limit Tc/W 0. 
It is reasonable to consider that NaxCo02-j/H20 is in the 
weak coupling region Tc/W <C since Tc ~ 5K. There- 
fore, we take the weak coupling limit from the beginning 
and solve the weak coupling model in eq. (3) within the 
mean field theory. 

In the following, we solve the mean field equation for 
Hi and find out the local minimum and saddle point of 
the free energy by performing the stability analysis. As 
a result, we search the minimum of the free energy and 
determine the phase transition. It should be noticed that 
the symmetry of Hamiltonian is reduced by the parallel 
magnetic field which we focus on. Then, we can only rely 
on the ^(1) gauge symmetry and there are 6x2—1 = 
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11 independent order parameters including the relative 
phase. 
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Fig. 1. Pairing state at H = and T = O.IT^''. Here, r° is the 
transition temperature at which the superconducting instabihty 
occurs, and we choose the parameter V so that = 6K. The 
notations of Pxy+, F-^y and Fz have been given in Table II in 
I.^*) In the Pxy+ + Fz-state, order parameters in the fxy+- and 
Fz-states coexist. 



Before showing the result under the magnetic field, we 
show the phase diagram at = and below in Fig. 1. 
We see the four states including the -Pxy+ + -FL-state. In 
this state, the p-wave order parameter with d- vector par- 
allel to the plane coexists with the /-wave order param- 
eter which has the d-vector along z-axis. This coexistent 
state is stabilized around the phase boundary between 
the p-wave and /-wave states so as to gain the conden- 
sation energy by making the superconducting gap more 
isotropic. In this paper, we denote this pairing state as 
p 4- /-wave state, for simplicity. In a small part of the 
parameter range, another pairing state such as Pz -f -Fxy- 
state may be stable. However, we ignore this possibility 
since the parameter region will be very small. 

As will be shown in the following sections, the pairing 
state under the magnetic field is determined by three ef- 
fects. (I) The spin-orbit coupling term which favors the 
pairing state shown in Fig. 1. (II) The magnetic field 
which is unfavorable for the c?- vector parallel to the field. 
(Ill) The momentum dependence of the superconducting 
gap, i.e. the more isotropic superconducting gap is, the 
more condensation energy is gained below T^. The mul- 
tiple phase diagram appears in the H-T plane as a result 
of these competing effects. 

3. Phase Diagram under the Magnetic Field 

In the following, we discuss the H-T phase diagram in 
each case, i.e. the p-wave, /-wave and p-l- /-wave states. 
Hereafter, we choose the parameter IJ so that — 6K 
consistent with experimental value, where T'^ is the tran- 
sition temperature of superconductivity at iJ = 0. Note 
that T;? is the maximum value of transition temperatures 
for Pxy+-, -Fxy- and Pz-states because the other state is 
not stabilized at T = Ti" and FT = in our calculation. ^ 



3.1 P-wave state 

First, we discuss the p-wave superconductivity which 
is realized in the large Jh region. We fix the parame- 
ters of two-orbital Hubbard model as Jb_/U = 7/24 and 
rie = 0.21 where we obtain T° = 6K by choosing U — 
6.47. In the linearized Dyson-Gorkov equation for the 
equal spin pairing, (eq. (1)), we obtain the wave functions 
0ix(A;) and 4>iyik) as A||(fc,i7rr) = -0ix(fc) -I- i(/)iy (fe) 
where A||(fc) is the eigenfunction corresponding to the 
maximum eigenvalue Aj^. In cq. (2), we obtain (j)2x{k) 
and (j)2y{k) from doubly degenerate eigenfunctions as 
A|x(fc,i7rr) = 02x(A;) and A|x(fc,i7rr) = 02y(fc). The 
difference between (j)i^{k) and 02x(fc) and that between 
0iy(fe) and 02y(fc) are slight because A ^ W. Because 
the amplitude of wave functons (j)ia{k) is not determined 
by the linearlized Dyson-Gorkov equation, we impose the 
normalization condition for (j)ia{k) as. 



ir.VM^^tanh^ 

2TP 



V (fc) 
k 



1, 2 and a = x,y), (4) 



where K is an arbitrary value. Although the coupling 
constant gi depends on the value of K , the physical quan- 
tities are obviously invariant for the choice of K . 

According to the symmetry of triangular lattice, we 
obtain gi = §2 and 174 = 35. Then, we have three inde- 
pendent parameters, namely f/i, g^ and 55. At A = 0, we 
obtain gi — 55 and 173 = 0. Thus, the role of spin-orbit 
coupling is represented by the two parameters, gi — 55 
and 33. 

The parameters gi -\- 53 and g^ are determined so as 
to reproduce the transition temperatures obtained in the 
linearized Dyson-Gorkov equation. From eq. (3), we ob- 
tain the gap equations at iJ = for the -Pxy+- and P^- 
states as. 



.,3)VM;)i^tanh.^^(^) 



k 



2T, 



(xy^ 



l=,aVM;^tanh.^^(^) 



k 



2T, 



(z) 



(5) 



(6) 



Here, Tc^^'^'' is the transition temperature for the Pxy+- 
state, which is determined by the criterion Aj^ = 1 at 
T = Tc^^^\ while Ti^^ is the transition temperature for 
the Pz-state and determined by the criterion A]-^ = 1 at 
T = T^^\ Since T^""^^^ > T^^^ at = as shown in 



T^^\ Since Ti"''+) > Pi^^ at H 

^(xy+. 



Fig. 1, we obtain Tr^+' = T° = 6K. We have obtained 
Pi'^^ = 0.865PO for A = 0.17 by solving eq. (2). Thus, 
91 + 53 and (75 are determined by the microscopic param- 



eters, such as U, Jh, and A, through Tc^^^^ and Pc 
respectively. 

In principle, the parameter 171 — 173 can be determined 
by the transition temperature for the Pxy--state, Tc^^~^ 
as. 



1. 2 M..„t.^>W 



k 



2E2{k) 2P,("^"^' 



(7) 



However, it is difficult to obtain the precise value of 
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Ti"'^ ^ from eq. (1) since ri""^ ^ 



< T^''y+\ Therefore, 



there remams an arbitrary parameter among gi and g^. 
This difficulty can be resolved by the microscopic analy- 
sis on the effective interactions V||(fc,A;') and F||(fc, /c'). 
As a result of the perturbation expansion in A, we have 
found that \gi — gz\ is much smaller than Igal.^^-* The 
former is the second order term in A, namely <7i — 35 = 
O(A^), while the latter is the first order term, namely 
(73 — 0(A). Therefore, we assume 51 — 55 = for a while 
and we will show the role of finite value of gi — 35 later. 




T/T 




TIT 



Fig. 2. The phase diagram in the p-wave superconducting state 
under the magnetic field along a;-axis. We choose A = 0.05 in (a) 
and A = 0.17 in (b). The notations of A, B, C and D phases are 
given in Table I. The solid (dashed) line shows the first (second) 
order phase transition. 



Figure 2 shows the phase diagram in H-T plane for 
two values of A. It should be noted that the effect of 
magnetic field is scaled by the typical value which is 
defined as ^isHp — ksT^- When we assume T° = 6K, we 
obtain Hp ^ 9T. Note that Hp is a slightly smaller value 
than the Pauli paramagnetic limit. We sec in Fig. 2 the 



four phases A, B, C and D whose notations are given in 
Table I. The coefficients a and (3 vary according to the 
magnetic field and temperature. We have not shown the 
small components in the order of 0{T^ /W) which are 
induced by the magnetic field and particle-hole asym- 
metry. Strictly speaking, A-, B-, and D-phase are the 
non-unitary state owing to these small components. Re- 
cently, it has been pointed out that the non-unitarity 
induced by the paramagnetic effect plays an important 
role for Sr2Ru04.^''^-' However, this effect is negligible in 
what follows since the spin-orbit coupling is much more 
significant for Na^CoOa • t/HaO than Sr2Ru04.^'^''*^^ 



A 


ap^x + pyy 


B 


ap^x + pyH ± Pp^z 


C 


(px±ipy)(ai Tiy) 


D 


(Px ± iPy)2 


E 




F 





Table I. The pairing states under the magnetic field. The first 
column shows the notation in this paper. The second column 
shows the d-vector. Note that 0<a, /3<1, ■|§^<0 and ^ > 0. 



The A-phase is connected to the Pxy+-statc ai H — Q 
in an adiabatic way. This phase is favored by the spin- 
orbit coupling and stabilized at low magnetic fields. The 
x-component of d-vector is reduced in the magnetic field 
so as to gain the Zeeman coupling energy, and therefore 
^ < 0. Note that the A-phase has a 2-fold degener- 
acy at = or at T = Tc{H) where d = ap^x + Pyi) 
is degenerate with d = apyX — p^y- This degeneracy is 
lifted for iJ ^ and T < Tc, and the former is stabilized 
because the order parameter in the py-wave symmetry 
4>iy{k) is more isotropic than that in the px-wave symme- 
try (j)ix{k). This hfting of degeneracy is very small and 
these two states are nearly degenerate. Therefore, it is 
expected that the orbital effect plays a more important 
role for the lifting of degeneracy. As will be discussed 
in §4.1, the degeneracy between d = ap^x + Pyi) and 
d = apyX — py^y is lifted by the the orbital effect, and 
the former is stabilized furthermore. According to the 
analysis on the Ginzburg-Landau free energy,"*^' "^^'^^^ the 
orbital effect on the lifting of this degeneracy is in the 
second order of order parameter, while the paramagnetic 
effect is beyond the forth order. 

The B-phase is stabilized at high magnetic fields. The 
reflection symmetry with respect to the a;j/-plane is vi- 
olated in the B-phase. In comparison to the A-phase, 
the Zeeman coupling energy is gained in the B-phase be- 
cause the x-component of d-vector is reduced, although 
the spin-orbit coupling term disfavors the B-phase. The 
A-B phase transition is described by the rotation of d- 
vector having the Px-orbital component, which is shown 
in Fig. 3 schematically. It should be stressed that the 
x-component of d-vector does not vanish even in the B 
phase owing to the g^-term in eq. (3). Because the am- 
plitude of Px-wave order parameter increases as increas- 
ing the magnetic field across the A-B phase transition, 
the superconducting gap becomes more isotropic in the 
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B-phase. In comparison to the C-phase, the B-phase is 
stabihzed at low temperatures since the superconducting 
gap is more isotropic in the spin space rather than the 
non-unitary C-phase. At high fields, the B'-phase where 
d = apyX — p^y ± PpyZ is nearly degenerate with the B- 
phase and can be stabilized. However, we have ignored 
this possibility since the B'-phase is disfavored by the 
orbital effect which is discussed in S4.1. 



A-phase 



B-phase 




Fig. 3. The schematic figure for the A-B transition. The d-vector 
with px-wave symmetry rotates towards z-axis. 



The C-phase is a non- unitary state, i.e. d x d* ^ 
where the time-reversal symmetry is broken. The C- 
phase can be stabihzed just below T — Tc{H) be- 
cause the d-vector is described by a linear combina- 
tion of two- fold degenerate state at T = Tc{H) as 
d — {apxX + PyV) ± i{apyX — p^y). In general, a linear 
combination of the degenerate states is stabilized just be- 
low Tc{H) so as to make the superconducting gap most 
isotropic. The superconducting gap is more isotropic in 
the spin space in the unitary A-phase, while that is more 
isotropic in the momentum space in the C-phase owing to 
the factor p^ ± ipy. As increasing the magnetic field and 
therefore decreasing a, the superconducting gap in the 
A-phase becomes anisotropic while that in the C-phase 
becomes isotropic. Therefore, the C-phase is stabilized at 
high fields. As decreasing the temperature, the B-phase 
is stable in comparison to the C-phase owing to the non- 
unitarity in the C-phase. Thus, the C-phase is stabilized 
at high magnetic fields and around T = Tc{H). 

The D-phase is the chiral superconducting state where 
the time-reversal-symmetry is broken. It has been ex- 
pected that this phase is stabilized in Sr2Ru04 at 
H = 0.^' However, the D-phase is not stabilized in 
NaxCo02 • J/H2O unless we assume a very small spin- 
orbit couphng as A = 0.05 in Fig. 2(a). 

In the following part, we clarify the physical proper- 
ties of each phase in order to compare with experimen- 
tal results and to suggest future experiments. We first 
show the temperature dependence of magnetic suscep- 
tibility X — M/H in Fig. 4. It is clearly shown that 
the magnetic susceptibility is smaller than the normal 
state value in both A-phase and B-phase owing to the 
a;-component of d-vector. The decrease of magnetic sus- 
ceptibility is reduced by increasing the magnetic field 
since the x-component of d- vector decreases so as to gain 
the Zeeman coupling energy. These behaviors are qualita- 
tively consistent with the recent NMR measurements in 
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Fig. 4. The temperature dependence of magnetic susceptibility, 
X = M/H. (a) A = 0.05 and H = O.lHp, 0.3Hp, 0.35Hp, O.bHp 
and Hp. (b) A = 0.17 and H = O.lHp, 0.5Hp, O.YHp, Hp and 
1.6Hp. 



the 0-site."' The magnetic susceptibility is almost tem- 
perature independent in the C-phase, but it decreases 
below Tr when we consider the orbital effect as in S4.1. 




- - X=0.05 

X=0.17 

■ - X=0.25 



H/H 



1.5 



Fig. 5. The magnetic field dependence of magnetic susceptibility 
at T = 0.1T°. We show the results for A = 0.05, A = 0.17 and 
A = 0.25. 



Figure 5 shows the magnetic field dependence of mag- 
netic susceptibility. As shown by the arrow, the mag- 
netic susceptibility rapidly increases through the A-B 
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phase transition. This kink will bo visible even if we con- 
sider the vortex state as in §4.1. Therefore, the detailed 
measurement on the magnetic field dependence of NMR 
Knight shift will be an interesting experiment. If wc as- 
sume A = 0.17 and T° = 6K as in Fig. 2(b), the A-B 
transition occurs around H = 0.7Hp ~ 6T. This mag- 
netic field is experimentally accessible. 



not incompatible with the existence of multiple phase 
transition. It is worth noting that the specific heat mea- 
surement has not detected the multiple phase transition 
in UPta at high fields. If we choose the magnetic field 
so that the second transition occurs near Tc{H) as in 
Fig. 6(b), the discontinuity in 7 becomes large. It may 
be possible to observe the A-B transition in the specific 
heat measurement by adjusting the magnetic field. 




ak x+k y+Bk z 



T/T' 



0.25 



5 0.15 




H/H 




— ak^x+kyV+pk^z 
.-. ak^x+kyV 



0.2 0.4 0.6 



T/T^ 

c 

Fig. 6. The temperature dependence of specific heat coefficient 
7 = C'/T (solid line). The dashed line shows the results when we 

assume the A-phasc in whole temperature range. We choose the 
magnetic field as (a) H = O.THp and (b) H = Hp. 



Since the A-B transition is the second order phase 
transition, the specific heat shows a discontinuity. We 
calculate the specific heat over temperature 7 = C/T. 
The typical behavior is shown in Fig. 6(a) where the 
discontinuity in 7 is not so significant. This is because 
the entropy is not so different between the A-phase and 
B-phase. It should be noticed that the A-B phase tran- 
sition is described by the rotation of d-vector as shown 
in Fig. 3 and the amplitude of order parameter changes 
only slightly. We note that the broadening of the spe- 
cific heat occurs in the vortex state. Therefore, it may 
bo difficult to observe such a small discontinuity in 7 
experimentally. Indeed, no evidence has been obtained 
for the multiple phase transition from the specific heat 
measurement in NaxCo02 • yH20.^''~^^) However, this is 



Fig. 7. The magnetization along z-axis over the magnetic field 
along i-axis. The parameters are chosen to be T = 0.1T° and 
A = 0.17. 



We find another interesting property of the A-B phase 
transition, namely the spontaneous rotation of the prin- 
cipal axis of magnetic susceptibility. In the A-phase, the 
principal axis is the same as the symmetry axis of crys- 
tal. On the other hand, the principal axis in the B-phase 
is misaligned from the crystal axis because the d-vector 
having the px-symmctry lies between the x and z direc- 
tion. This particular property is illuminated by showing 
the magnetization along the z-axis, under the paral- 
lel magnetic field. Although the magnetization is zero 
in the A-phase and in the normal state, M-^ has a finite 
value in the B-phase due to the off-diagonal element of 
susceptibility tensor. We show the magnetic field depen- 
dence of \M.^\/H^ in Fig. 7. The absolute value is shown 
because the sign of is different between the 2-fold de- 
generate states in B-phase (see Table I). It is shown that 
the magnetization \M.^\/H^ shows a sharp peak in the 
B-phase close to the A-B transition. When we increase 
the magnetic field from the A-B transition line, |Mz|/i?x 
first increases owing to the increase of f3 and then de- 
creases owing to the decrease of a. Since the appearance 
of |Mz| is a particular property of B-phase, the measure- 
ment of this quantity will be an interesting experiment. 
We note that \Mz\ is finite also in the C-phase because 
of the non-unitarity. However, the value in the C-phase 
is in the order of /W and therefore much smaller than 
that in the B-phase. 

Thus far, we have assumed 31—55 = 0. At the last 
of this subsection, we discuss the role of the spin-orbit 
coupling through the finite value of gi — 5(5. We fix A = 
0.17 as Fig. 2(b) and choose two values of gi —55 so that 
Tc2 = 0.825T° and = 0.905T°. Here, Tc2 is defined 
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Fig. 8. The phase diagram when the value of g\ — gs is finite. We 
fix A = 0.17 and choose T^i = 0.825T0 (31 - < 0) in (a) and 
Tc2 = O.gOSrO (gi - 95 > 0) in (b). 



as, 



k 



2E2{k) 



2T, 



c2 



(8) 



Since = ri'^ = QMhT^ for gi -55 = 0, gi -35 < in 
the former case and gi—g^ > in the latter case. Figure 8 
shows the phase diagram in each case. In both cases, 
the phase diagram in the low field region is similar to 
Fig. 2(b). For instance, the A-B phase transition occurs 
around H — Q.7Hp. On the other hand, the pairing state 
at high magnetic fields is affected by the finite value of 
gi — 55. The D-phase is stabilized when gi — gz < 
while the C-phase is stabilized when gi — gs > 0. These 
results are understood by noticing that ga-term in eq. (3) 
is ineffective at high magnetic fields since a 1. Then, 
the D-phase is favored when the (75-term is large while 
the C-phase is favored when the gi-term is large. 

As shown in 1,^^^ the second order term in A stabilizes 
the Pxy-i--state furthermore. This means that gi — g^ > 



0. Since a large value of gi — 55 has been assumed in 
Fig. 8(b), the role of gi — g^ may be overestimated. We 
expect that the C-phase is slightly favored by the finite 
value of gi — 35 in NaxCo02 • 2/H20. 




H/H 



1.5 



Fig. 9. The magnetic field dependence of magnetic susceptibility 
at T = O.ITO and A = 0.17. We show the results for gi - gs < 0, 
91 ~ 95 = S'lid 91 ~ 95 > Oj respectively. 



Figure 9 shows the field dependence of magnetic sus- 
ceptibility for gi— g5 ^ 0. Although the behaviors in the 
low field region are independent of the value of gi— 55, the 
magnetic susceptibility at high fields is affected. The dis- 
continuity appears at the B-D and B-C transition which 
are the first order phase transition. The magnetic suscep- 
tibility in the C- and D-phases does not decrease from 
the normal state value. We note again that the magnetic 
susceptibility decreases in the C-phase when we take into 
account the orbital effect as in §4.1. 

3.2 F-wave state 

The phase diagram in the /-wave superconducting 
state is much simpler than the p-wave superconducting 
state. The multiple phase transition is not expected un- 
der the magnetic field along the plane even if the Fy^y- 
state is stabilized at H = 0. Among the 2-fold degeneracy 
in _Fxy-state, d = af2X + fiy is favored by the magnetic 
field along x-axis, as shown in Fig. 10 schematically. Al- 
though the magnetic susceptibility slightly decreases ow- 
ing to the ^-component of d- vector, the decrease is much 
less than 1% in our estimation because a <^ 1.^'*-' Then, 
the level crossing between the Fxy- and i^z-state due to 
the Zeeman coupling energy occurs at the extraordinary 
high magnetic field which will be higher than Hc2- 



H=0 



H//X 



f^x-afgy 
afsX+f^y 



Fig. 10. The schematic figure for the energy level in the /-wave 
superconducting state. 



If we apply the magnetic field along z-axis, the rota- 
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tion of c?-vector occurs when the Fz-state is stabihzed at 
H — 0. However, it will be difficult to observe this phase 
transition because the magnetic field is very small owing 
to the remarkably small splitting of Tc between i^xy- and 
Fz-states, as shown in I.^^^ 

From these discussions, it is expected that the single 
phase would be observed in the /-wave superconducting 
state in both directions of magnetic field. We note that 
almost temperature independent magnetic susceptibility 
expected in the /-wave state is incompatible with the 
resuhs of Co-NMR9'1°'5i,52) ^^^^ O-NMR.i^' 

3.3 P + F-wave state 

Finally, we investigate the coexistent p + /-wave state. 
We denote this state as E-phase in Table I. In the Jh- 
Uc phase diagram in Fig. 1, the coexistent state is sta- 
bilized in the region denoted as Pxy+ + F^. Although 
such a coexistent state with different orbital symmetry 
is usually realized only accidentally, the Pxy+ + -F^-state 
is stabilized in a relatively large parameter region. This 
is because the p-wave and /-wave superconductivities are 
nearly degenerate. We have provided a clear explanation 
on this degeneracy by analyzing the orbital character of 
six hole pockets.*^-' 

In order to investigate the p + f-wave state in the weak 
coupling approximation, we consider the gi-, 52- and (73- 
terms having the p-wave symmetry and (75-term having 
the /i-wave symmetry. Hence, gi — 32 due to the sym- 
metry of triangular lattice, however g^ ^ 175 because the 
/i-wave symmetry is not degenerate with the the /2-wave 
symmetry which is induced the g4-term. Hereafter, we ig- 
nore the (74-term since the for the /2-wave symmetry is 
very small. For simplicity, we assume the orbital wave 
functions in the p^,-, py- and /i-symmetry as 0ix(fc) = 
■\/3sin ^fcx cos ^fcy, (/'iy(fc) = sinfcy -I- sin ifcy cos ^fcx 
and (t>2y{k) = sin ifcy (cos ifcy — cos^fcx), respectively. 
These wave functions denote the Cooper pairing between 
the nearest neighbor sites. The parameters gi — g2, 
gs are determined from eqs. (5) and (8) by assuming 



tc2 



0.865Ti''''+' and ri''''+' = 6K. We choose the crit- 
ical temperature in the /i-wave symmetry, namely Tc 
as a parameter and determine 55 by eq. (6), although 
Ti^' can be determined by the microscopic parameters, 
such as Jh and rio- 

We find three kinds of the H-T phase diagram which 
are determined by Tc^\ These are classified as (A) 



(z) 



< Te2 < Ti"^+\ (B) Te2 < T^' < T 
'(xy+) , rp(z) 



.(z) 



(xy-H) 



and 



(C) Tc2 < ' < n '■ Since the case (B) is obtained 

in the largest parameter range of two-orbital Hubbard 
model, we assume the case (B). The phase diagram in 
this case is shown in Fig. 11. 

If we apply the magnetic field along x-axis, there ap- 
pears a tetracritical point as in Fig. 11(a). The A-phase 
is stabilized around T ~ Tc at low magnetic fields since 
Ti^^ < Ti''y+^ At high magnetic fields, the F-phase (F^- 

state) is stabilized round T = TdH) since Tc2 < Ti^\ 
In the low temperature region, the E-phase is stabilized 
in comparison to the A-phase and F-phase since the su- 
perconducting gap is isotropic. The coexistent E-phase 
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Fig. 11. The H-T phase diagram in the p + /-wave supercon- 
ducting state for (a) || a; and (b) H \\ z. Wo choose the pa- 
rameters, gi = §2, 93 and 35, so that Tc2 
T<?"> =0.95Ti"^+'. 



0.865Ti''^+' and 



appears below the second transition temperature which 
is lower than Tc . 

We note that Fig. 11(a) is obtained by neglecting the 
possibility of Fxy-state. The Fxy-state can be stabilized 
instead of the F-phase when the Fxy-state has higher 
transition temperature than the Fz-state. In this case, it 
is sufficient to substitute the Fxy-state for the Fz-state in 
the following discussions. Note that the B-phase and D- 
phase are not stabilized in the case (B) since Tc2 < Tc^\ 

Fig. 11(b) shows the phase diagram under the mag- 
netic field along 2-axis. In this case, we obtain a = 1 
in Table I since the symmetry of point group is not 
reduced. It is shown that the E-phase disappears at 
high magnetic fields because the magnetic field disfa- 
vors the d-vector along z-axis. The C-phase is stabilized 
around T = Tc{H). The degeneracy in the C-phase is 
lifted because the magnetization arising from the non- 
unitarity couples to the magnetic field. In the present 
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case, d = (px+iPy)(i — iy) is favored, where only the down 
spins are superconducting. This is because the DOS at 
Fermi level is larger for down spins owing to the van-Hove 
singularity at the Brillouin zone boundary. 




I M I I I I I I I I I L 

• 0.2 0.4 0.6 0.8 1 



T/T° 

c 

Fig. 12. The temperature dependence of magnetic susceptibility 
(a) along i-axis and (b) along z-axis. We choose the same pa- 
rameters as in Fig. 11 and the values of magnetic field are shown 
in the legend. The arrows show the second phase transition. Tae 

is the transition temperature from A- to E-phase at H = QAHp 
and Tbf is that from E- to F-phase at i/ = H^. 



One of the characteristic properties in E-phase is the 
decrease of magnetic susceptibility in all directions of 
magnetic field. This is because the rf-vector has all com- 
ponents, i, y and z. We show the temperature depen- 
dence of magnetic susceptibility along a^axis in Fig. 12(a) 
and that along z-axis in Fig. 12(b). This feature of E- 
phase is interesting because the decrease of NMR Knight 
shift along a;-axis has been reported by Kobayashi et 
al}'^^ This consistency should be contrasted to the p- 
wave state (§3.1) and /-wave state (§3.2) where the mag- 
netic susceptibility along Zraxis does not decrease below 
Tc. We note that the magnetic susceptibility does not 
decrease in the A- and C- phases in Fig. 11(b). 

Another important feature of the p+f coexistent state 
is the existence of double transition at zero magnetic 
field. This is in sharp contrast to the phase diagram in 
the p-wave state (Fig. 2) which shows no double transi- 
tion at if = 0. The second transition from A to E-phase 
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Fig. 13. The temperature dependence of the specific heat coeffi- 
cient in the p + /-wave state for two values of T^f^ . The solid 
and dashed lines show the results for T^:''^ = 0.95Tc^^^' and 
ri^' = 0.9Ti''^+\ respectively. We fix T^2 = 0.865ri''^+' and 
i? = 0. 



induces a clear discontinuity of specific heat coefficient 7 
as shown in Fig. 13. This is because the superconduct- 
ing gap is much more isotropic in the E-phase than in 
the A-phase. The p-wave order parameter is remarkably 
anisotropic on the Cg-Fcrmi surface since both p^- and py- 
wave order parameters vanish around the K-point owing 
to the periodicity of Burillouin zone. In contrast to that, 
the /i-wavc order parameter is almost isotropic on the 
Cg-Fermi surface. Therefore, the entropy is significantly 
reduced by the appearance of /-wave order parameter. 
Although the specific heat has been measured at zero 
magnetic field, no sign of the second phase transi- 
tion has been reported up to now. We suggest that an 
experimental search for this second phase transition will 
be an interesting future issue. 

We have calculated the NMR 1/TiT in the p /-wave 
state, but we can sec a very weak anomaly at the second 
phase transition. Therefore, it will be difficult to ob- 
serve this second phase transition by the measurement 
of NMR 1/TiT. 

4. Roles of Vortex State and Aig-Orbital 

In §3, we have ignored two effects which may affect 
the H-T phase diagram. One is the orbital effect which 
induces the vortex state. The other is the role of aig- 
orbital. In this section, we investigate these effects within 
the qualitative discussion. 

Vortex state 

First, we discuss the vortex state. Most of the uncon- 
ventional superconductors including NaxCo02 -2/1120 are 
classified into type II where the vortex state is stabilized 
at Hci < H < Hc2- The spatial dependence forming 
the vortex state is not important for the d-vector, qual- 
itatively when the orbital part of order parameter has 
one component. The /-wave state discussed in §3.2 is 
the case. However, the rf-vector can be affected by the 
orbital effect when the orbital part has multi-component 
as in the p-wave and p-|-/-wave states. In these states, the 
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d-vector in the vortex state is discussed in the following 
way. 

We discuss the vortex state on the basis of the 
Ginzburg-Landau theory which provides qualitative and 
clear understandings. In this theory, the free energy is 
given by the functional of order parameters as, 

F = Fo + Fg, (9) 

where Fq is the uniform term and Fq is the gradient 
term. Here, we define the order parameters d'^^ through 
the (i-vcctor as d = (i2'Vi/3(fc)Q: where i = 1 for 

a = X, y and i = 2 for a = z. In this subsection, we 
analyze the GL functional within the quadratic order and 
determine the pairing state at T = Tc{H) and discuss the 
d- vector below Tc{H). 

We first consider the p-wave superconductivity. Within 
the quadratic order of uniform term Fq , there is a 2-fold 

degeneracy between (di^^d^""') ^i(7lfs5' "Tifc^) 
and (di"\d(^)) = A^i^^,^^). These states are 
obtained from the doubly degenerate Pxy+-state in an 
adiabatic way and nearly degenerate in the A-phasc. By 
taking into account the two component order parame- 
ters, namely Ai and A2, the quadratic order term in Fq 
is written as Fq = a{T - T^""^^ {H))[\Ai\'^ + IA2P] where 
Tc^^\h) is the transition temperature obtained in §3.1. 
Note that the A-phasc and C-phase are described by a 
linear combination of these two component order param- 
eters. We obtain (A1.A2) cx (0.1) in the A-phase, and 
(Ai, A2) cx (±i, 1) in the C-phasc. 

When we consider the magnetic field along a^axis and 
ignore the spatial dependence in a^direction, the gradient 
term is given as, 

Fg = kI{\dA^^\^ + |i?yd(-)|n + KtlDydi^^^ 
+K^^\{\Dydiy^' + |Dyd(^)p} + K^\DySy^\^ 
+i^|{|D,d(-)|2 + |D,d(-)|2 + |D,d(^)|2 + |D,d(^)|2} 

+K^{\DA^^\' + \DA'^f}, (10) 

where Dj = dj — 2\eAj is the gauge invariant differ- 
ential operator. The coupling constants are obtained in 
the clean limit as, k\ : kI : kI = < \4>i^(k)\'^v'^ > : 
< |</)iy(fc)P^;2 > : < |0i,(fc)|2i;2 > and : : K^r 
= < |02.(fe)|2,;2 > : < |^2y(fe)P«2 > : < |</,2.(fe)|2«2 > 
where Va is the Fermi velocity along a-direction and <> 
is the average on the Fermi surface. We obtain Kf ^ 
in general, but these are nearly the same values because 
A <C W. Therefore, we simply assume Kf = = Ki. 
If the cylindrical Fermi surface and purely p-wave order 
parameter are assumed, wc obtain K2 — iKi. However, 
we find that K2 < Ki for the Fermi surface relevant for 
NaxCo02 • yH20. 

It is simply understood that the 2-fold degeneracy be- 
tween Ai and A2 is lifted by taking the gradient term. 
By taking into account the spatial dependence of order 
parameters as Ai Ai(r) and A2 = A2(r) and solving 
the linearized GL equation, we obtain two different Tc 
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r(xyl)(g) ^ r(xy)(^) _ |^ "'f; + f ^ ^g, (11) 

T(xy2) ^H) = T(-) (H) - ^ J^^^^sH. (12) 

ay 1 + 

The former is the critical temperature for the order 
parameter Ai, while the latter is that for A2. Since 
K2 < Ki and a < 1, wc obtain 

T.l'^yi) {H) < T^""^^^ {H). (13) 

Thus, the degeneracy between d = ap^x + pyfj and d = 
apyX — p^y at T = Tc{H) is lifted by the orbital effect 
and the former is stabilized. Therefore, the A-phasc is 
always favored rather than the C-phase at T = Tc{H) 
in contrast to Fig. 2(b). The C-phase can be stabilized 
below Tc{H) through the second order phase transition. 
This is the most significant role of the orbital effect for 
the d-vector. 

As for the d-vector along z-axis, the degeneracy be- 
tween d = pxZ and d ^ pyZ is lifted by the gradient term 
in the same way. We obtain the Tc as, 

Ti^i) {H) = Ti^) (H) - ^ ^/kJ{^H, (14) 
Ti^'^) (H) = Ti^) (H) - ^ ^/K^H, (15) 

where Tc^^ {H) is the transition temperature of D-phase 
without including the orbital effect. We obtain 

T^^^\H) < T^^^\H), (16) 

and therefore d = pyZ is more favorable than d = p^z. 
This splitting of transition temperature is due to the dif- 
ference of coherence length along y-axis between the Px- 
andpy-states. The same mechanism for the splitting of T^. 
has been discussed for Sr2Ru04.'''^^ Then, Agterberg has 
predicted the double transition a.t H > iJ^i- However, 
the experimental observation in Sr2Ru04 is incompati- 
ble with the theoretical prediction. Recently, Udagawa 
has proposed another idea which resolves this inconsis- 
tency.^^) 

Whether T^^^\h) < T^"'^^\h) or not depends on the 
parameters such as Ki and T^''^\h) - T^^\h). When 
we assume gi — 95 = as in Fig. 2, we always obtain 
T^^'^\h) < T!r''^\H). Then, the A-phase is always sta- 
bilized at T = T^{H). 

We show the schematic figure of H-T phase diagram in 
Fig. 14(a-c). The C-phase appears below Tc{H) when we 
assume Hc2 3> Hp. Then, the phase diagram is expected 
as Fig. 14(a). When Hc2 Hp, the C-phase may disap- 
pear as in Fig. 14(b). The B-phase is also destabilized 
by the gradient term at high magnetic fields, however 
stabilized at low magnetic fields. When Hc2 is decreased 
furthermore, the B-phase disappears as in Fig. 14(c). In 
this case, another phase d = pyZ can be stabilized at 
high magnetic fields if we assume (?i — (75 < 0. However, 
this is not expected in our results for the two-orbital 
Hubbard model. ^"'^ Unfortunately, there are two differ- 
ent experiment reports for the i?c2 in NaxCo02 • yH20. 
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(a) 

H 




(b) (c) 




Fig. 14. The schematic figure of H-T phase diagram in the p- 
wave superconducting state when we take into account the vortex 
state. The magnetic field is applied along s-axis. The parameter 
Hc2/Hp is decreased from (a) to (c). 
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from the 2-fold degeneracy in the Pxy+-state which is re- 
lated to the p-wave state (§3.1) and p + /-wave state 
(§3.3). 

The degeneracy in the Pxy+-state is due to the C/(l) 
symmetry corresponding to the in-plane rotation of 
spin.^^-* This symmetry is broken in the three-orbital 
Hubbard model including the aig-orbital because the in- 
plane component of orbital moment has the matrix ele- 
ment. Therefore, the degeneracy in the two-orbital Hub- 
bard model will be lifted at zero magnetic field. However, 
it is expected that the lifting of degeneracy is small be- 
cause the flig-orbital is not so important for the super- 
conducting instability.'^''''*^' Then, the phase transition 
can occur in the low magnetic field region. 




J 



A' «PyX-Pxy'. 
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The resistivity measurement has reported Hc2 compara- 
ble to Hp,^°^ while the magnetization measurement has 
reported Hc2 much larger than Hp.^^^ 

The orbital effect in the p + /-wave state is qualita- 
tively different from the p-wave state in the sense that 
Kf 7^ K;^. In this case, there are two possibilities of the 

high field phase diagram. When " f+a'^^ ^3 > -^2~^3~' 
the phase diagram is qualitatively the same as Fig. 11(a). 
Since the gradient term favors the /-wave order param- 
eter, the tetracritical point shifts to the lower magnetic 

field. If we assume " f+a'^^ < ^i'K'tj the /-wave or- 
der parameter is disfavored by the gradient term. Then, 
the E-phase disappears in the high field region, and 
the phase diagram is similar to Fig. 11(b). The micro- 
scopic examination of the magnitude relation between 

" i^a-^^ ^3 K^K^ is difficult because it is neces- 
sary to estimate ifj and Kj^ originating from the dis- 
persion relation along z-axis. 

In this paper, we have focused on the phase transition 
related to the d-vector and we have not discussed the 
phase transition in the vortex configuration. Because the 
vortex lattice structure is sensitive to the pairing state, 
the latter phase transition may occur in addition. We 
have simply ignored this possibility because the mag- 
netic and thermodynamic properties are not sensitive to 
the vortex lattice structure. We think that the more so- 
phisticated and quantitative calculation on the vortex 
state is one of the fascinating future issues. 

4-2 Role of aig- orbital 

We suggest another possibility of phase transition in- 
duced by the aig-orbital. This phase transition originates 



Fig. 15. The schematic figure of H-T phase diagram when the ad- 
ditional phase transition is induced by the aig-orbital. Although 
we have assumed the p-wave superconducting state, qualitatively 
the same phase transition can occur in the p + /-wave state. 



Let us consider the magnetic field along a;-axis. Then, 
not only the paramagnetic effect (§3) but also the orbital 
effect (§4.1) favor d = ap^x + Pyjj, i-e. A-phase. There- 
fore, the level crossing occurs under the magnetic field 
when the other Pxy+-state, i.e. d = PyX — py,y, is stabi- 
lized at H — 0. We show the schematic phase diagram in 
Fig. 15 where the additional phase transition is assumed. 

The experimental identification for this phase tran- 
sition may be difficult because the magnetic and ther- 
modynamic properties are quite similar between these 
states. However, we note that this phase transition is as- 
sociated with the lattice distortion through the electron- 
phonon coupling. Therefore, this transition may be de- 
tected by the measurements such as the magnetostric- 
tion, thermal expansion and ultrasonic attenuation and 
so on. Actually, the multiple phase transition in UPta 
has been identified by these methods. ^°'^^' 

5. Summary and Discussion 

In this paper, we have investigated the multiple super- 
conducting phase transition in the two-orbital Hubbard 
model representing the Cg-orbitals in NaxCo02 • J/H2O. 
We have derived the weak coupling model from the lin- 
earized Dyson-Gorkov equation within the second order 
perturbation theory and solved it within the mean field 
theory. 

In order to determine the pairing state under the mag- 
netic field, we take into account the atomic spin-orbit 
coupling term as well as the Zeeman coupling term. As a 



J. Phys. Soc. Jpn. 

result of these competing effects, the multiple phase tran- 
sition occurs in the superconducting state. We have ex- 
amined the phase diagram in the p-wave superconductiv- 
ity, /-wave superconductivity and their coexistent state, 
respectively. The obtained phase diagrams are qualita- 
tively different because the role of spin-orbit coupling is 
quite different between the p-wave and /-wave order pa- 
rameters. The multiple phase transition occurs under the 
magnetic field along the plane in the p-wave supercon- 
ducting state, while that occurs at zero magnetic field in 
the p + f-wa.ye state. However, the multiple phase tran- 
sition does not occur in the /-wave state. The obtained 
results are summarized in Table II. 



Symmetry 


p-wave /-wave p + /-wave 


Double transition 


H =/=0 none H = 


K{H ± z) 


O X o 


K{H II z) 


X X o 



Table II. Summary of the obtained results. The first column 
shows the pairing symmetry obtained in the multi-orbital Hub- 
bard model. The second column shows the possibility of multiple 
phase transition. The third and forth columns show the Knight 
shift along the plane and along the z-axis, respectively. The sym- 
bol O ( X ) means that the Knight shift decreases (does not 
decrease) below Tc. 



The characteristics of each pairing state have been il- 
luminated by showing the magnetic susceptibility. When 
we apply the magnetic field along x-axis, the magnetic 
susceptibility decreases in the p-wave state and in the 
P + /-wave state, while that does not decrease in the 
/-wave state. The decrease of in-plane Knight shift has 
been reported by the Co-NMR9'1°'5i,52) ^^^^^ O-NMR^^) 

which is consistent with the p-wave state and p + /-wave 
state. When we apply the magnetic field along z-axis, 
the magnetic susceptibility decreases in the p + f coex- 
istent state, while that does not decrease in the p-wave 
and /-wave states. A report of c-axis Knight shift in the 
Co-NMR^''^ is consistent with the former. 

Prom these observations, the p -I- /-wave state is ab- 
solutely consistent with the Knight shift measurements. 
However, it is still difficult to conclude the pairing state 
which is consistent with the experimental results in a 
comprehensive way. A difficulty of the p + /-wave state is 
that the second phase transition has not been observed at 
zero magnetic field. For instance, specific heat^^"^^^ mea- 
surements have not detected any anomaly. Another dif- 
ficulty is the power-law behaviors observed in the NMR 
1/TiT and specific heat. Since the superconducting gap 
is almost isotropic in the p + /-wave state, the exponen- 
tial behaviors will be observed in the low temperature 
region. 

Contrary to the p + /-wave state, the p-wave super- 
conducting state is consistent with many experimental 
results except for the c-axis Knight shift. For instance, 
the nearly power-law behaviors are expected because the 
superconducting gap is very anisotropic. This is partly 
because both px- and Py-wave order parameters van- 
ishes around the K-point near the Cg-Fermi surface, and 
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furthermore because the nodal quasi-particles exist in 
the aig-Fermi surface. The p-wave state is also consis- 
tent with the nSR measurement which does not de- 
tect any spontaneous time-reversal-symmetry breaking 
in contrast to Sr2Ru04.^^^ 

It should be noted that there remain some issues to be 
resolved when we assume the spin singlet pairing. As for 
the d-wave pairing, the /uSR measurement is incompati- 
ble with the dxy +idx2-y2-wave state which is expected in 
the triangular lattice. ^^"^^^ As for the i-wave pairing, 
it seems to be difficult to find a microscopic mechanism 
which induces such a high angular moincntum pairing 
with Tc ~ 5K. 

From these discussions, we think that further devel- 
opments in the theory and experiment are required to 
identify the pairing symmetry. The experimental search 
for the multiple phase diagram is particularly interesting 
because the phase transition related to the d-vector is a 
characteristic feature of spin triplet superconductivity. If 
the multiple phase transition is discovered in the super- 
conducting state, that will be a conclusive evidence for 
the spin triplet pairing. 

In order to suggest future experiments, we have clari- 
fied the nature of multiple phase transition. In the p-wave 
superconductivity, the A-B transition is promising to be 
observed. It is shown that the in-plane Knight shift shows 
a kink at the A-B transition. We find that the magneti- 
zation along the c-axis appears in the B-phase under the 
in-plane magnetic field. This spontaneous rotation of the 
principle axis of susceptibility tensor is a characteristic 
phenomenon in the p-wave superconducting state. 

In order to conclude the p -|- /-wave superconductivity, 
the second phase transition at zero magnetic field should 
be observed experimentally. The accurate measurement 
of specific heat around T =Tc is highly expected because 
the second transition can be masked by the broadening 
effect arising from the disorder. 

The measurements on the crystal lattice structure 
are also interesting since the multiple phase transi- 
tion is associated with the distortion of lattice through 
the electron-phonon coupling. For instance, the magne- 
tostriction, thermal expansion and ultrasonic attenua- 
tion have played an important role in the studies on the 
heavy fermion superconductors. 

In summary, we have examined the possibilities of mul- 
tiple phase diagram under the parallel magnetic field and 
suggested some future experiments. We hope that the 
pairing state in NaxCo02 • yH20 will be resolved by fu- 
ture experiments on the multiple phase transition. 
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